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We use magnetotunnelling spectroscopy as a non-invasive 
probe to produce two-dimensional spatial images of the proba- 
bility density of an electron confined in a self-assembled semi- 
conductor quantum dot. The images reveal clearly the ellipti- 
cal symmetry of the ground state and the characteristic lobes 
of the higher energy states. 



I. INTRODUCTION 

Quantum dots (QDs) are characterised by relatively 
small number of electrons confined within an island with 
a nanometer dimension. They can confine the motion 
of an electron in all three spatial dimensions The 
strong confinement in the QD gives rise to a set of dis- 
crete and narrow electronic energy levels similar to those 
in atomic physics. The epitaxial growth of lattice mis- 
matched InAs on GaAs or AlAs opens new possibilities 
for the simple fabrication of semiconductor nanostruc- 
tures. InAs QDs are formed in situ during growth due 
to the relaxation of a strained InAs wetting layer on 
GaAs or AlAs The particular interest lies in their 
uniformity and small size: a lateral dimension 10 — 20 nm 
and height 3 — 4 nm. Several theoretical approaches have 
been used to calculate the eigenstates of InAs QDs [||. 
The results of the calculation depend strongly on the 
assumed shape and composition of the QDs. Experimen- 
tally, the quantized energy levels of a given potential can 
be probed using various spectroscopic techniques. The 
corresponding wave functions are much more difficult 
to measure. Information about the extent of the car- 
rier wavefunction for the ground state of a QD was ob- 
tained from tunneling measurements in a magnetic field 
Q. Also, the anisotropy of electronic wave function in 
self-aligned InAs QDs was deduced from magnctic-ficld- 
dependent photoluminescence spectroscopy Q. How- 
ever, until recently there have been no reported measure- 
ments of the detailed spatial form of the wavefunctions of 
the ground and excited states of the QD. Recently, it has 
been demonstrated that magnetotunneling spectroscopy 
can be employed as nonivasive probe to produce images 
of the probability density of the electron confined in a QD 



. In this work, we use magnetotunneling spectroscopy 
(MTS) to investigate in detail of the spatial form of the 
wave function of the electron states of a double-barrier 
resonant tunnelling diode with InAs QDs embedded in 
the centre of a GaAs quantum well. We measure the 
dependence of the resonant tunnelling current through 
the QD states as a function of magnetic field, B, ap- 
plied perpendicular to the tunnelling direction. This al- 
lows us to map out the full spatial form of the probabil- 
ity density of the ground and excited states of the QDs. 
The electron wavefunctions have a biaxial symmetry in 
the growth plane, with axes corresponding quite closely 
(within measurement error of 15°) to the main crystal- 
lography directions X - [Oil] and Y - [233] for (311)B 
substrate orientation. For a similar InAs QD structure 
grown on a (lOO)-substrate we also obtained characteris- 
tic probability density maps of ground and exited states. 

II. EXPERIMENTAL DETAILS 

The InAs QDs are embedded in a n — i — n, resonant 
tunnelling diode. The samples were grown by molecu- 
lar beam epitaxy on a GaAs substrate with the orien- 
tations (100) and (311)B. A layer of InAs QDs, nom- 
inally 2.3 monolayer thick, was placed in the centre of 
a 9.6 nm wide GaAs quantum well (QW) with 8.3 nm 
Alo.4Gao.6As confining barriers, sandwiched between two 
nominally-undoped 50 nm GaAs spacer layers. The in- 
trinsic region is surrounded by graded n-type contact 
layers, with the doping concentration increasing from 
2 x 10 17 cm 3 to 3 x 10 18 cm -3 . The layers were grown 
at 600° C, and there was a growth interrupt before the 
QDs were grown at 480 °C. For comparison, we also stud- 
ied two control samples grown with the same sequence of 
layers, except one has only a thin InAs two-dimensional 
wetting layer (i.e. containing no QDs) and the other has 
no InAs layer at all. The samples were processed into 
circular mesa structures of diameters between 50 nm and 
200 fim, with ohmic contacts to the doped regions. 
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III. RESULTS AND DISCUSSION 
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FIG. 1. Schematic conduction band profile under an ap- 
plied bias of an n-i-n GaAs/(AlGa)As double barrier resonant 
tunneling diode incorporating InAs self-assembled quantum 
dots (QDs). Inset: Orientation of the magnetic field, B, and 
current,/, in the magnetotunnelling experiment. X and Y 
define the two main crystallographic axis, [011] and [233], re- 
spectively, in the (311)-oriented GaAs substrate, a and j3 
indicate, respectively, the direction of B and of the momen- 
tum acquired by the tunnelling electron due to the action of 
the Lorentz force. 

Figure 1 shows a schematic energy band diagram for 
our device under bias voltage. X and Y define the two 
main crystallographic axis in the plane perpendicular 
to growth direction, Z (see inset). The layer of InAs 
QDs introduces a set of discrete electronic states below 
the GaAs conduction band edge. At zero bias voltage, 
equilibrium is established by electrons diffusing from the 
doped GaAs layers and filling some dot states. The re- 
sulting negative charge in the QW produces depletion 
layers in the regions beyond the (AlGa)As barriers. By 
applying a bias voltage to the emitter layer, V, the QD 
energy level is shifted in energy with respect to both con- 
tacts. When a particular dot state is resonant with an 
adjacent filled state in the biased electron emitter layer, 
electrons tunnel through the dot into the collector and a 
current flows as shown schematically in Fig.l. Therefore, 
as we adjust the voltage we can study different energy 
states of the QDs. At sufficiently high voltages we are 
able to observe two separate resonances in the current 
related to confined subbands of the QW states. 
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FIG. 2. Low-temperature (T = 1.2 K) current-voltage 
characteristics, I(V). Dots are grown on (311)B substrate 
orientation. 

Figure 2 shows a typical low-temperature (T = 1.2 K) 
I(V) curve. The device contains InAs QDs grown on a 
(311)B-oriented GaAs substrate. Similar were obtained 
for QDs grown on a (lOO)-substrate. We observe a series 
of resonant peaks corresponding to carrier tunnelling into 
the dot states. Pronounced current features appear at 
biases as low bias as 55 mV. They are superimposed on a 
rising background current and cannot be resolved for V > 
200 mV. These features are not observed in our control 
samples which do not contain QDs and are thereforewe 
can ascribe them to the InAs QD layer. Despite the large 
number of quantum dots in our sample (10 6 — 10 for a 
100/jto diameter mesa), we observed only a small number 
of resonant peaks over the bias range (~ 200 mV) close 
to the threshold for current flow. This behaviour has 
been reported in earlier studies J|j6]]^-[l2j and, although 
not fully understood, is probably related to the limited 
number of conducting channels in the emitter that can 
transmit electrons from the doping layer to the quantum 
dots at low bias. There is no reason to believe that the 
dots studied are atypical of the distribution as a whole. 
On increasing the temperature to 4.2 K, the main peaks 
are still prominent, but much weaker features, which may 
be related to density-of-state fluctuation in the emitter 
Jl3| , are strongly suppressed. A key observation is that 
many peaks look similar so we cannot tell if the peaks 
are due to tunneling through the states of a single dot 
or several dots. In the following, we will concentrate on 
three voltage regions labelled ("A"), ("B") and ("C"). 
We will focus on the magnetic field dependence of the QD 
resonances and on how this provides detailed information 
about the form of the wavefunction associated with an 
electron in the QD state. 



2 




-8-6 -4 -2 2 4 6 8 -8 -6 -4 -2 2 4 6 8 



B (T) B (T) 

FIG. 3. (a) Low-temperature (T = 4.2 K) I(V) character- 
istics in the presence of a magnetic field, B. The direction of 
B is perpendicular to the current flow, (b) and (c) Depen- 
dence of conductance, G, on magnetic field for B parallel to 
[Oil] for different QDs states. B was varied from to 8 T 
with the step of 0.5 T. 

Figure 3(a) shows the low-temperature (T = 4.2 K) 
I(V) characteristics in the presence of a magnetic field, 
B. The direction of B is perpendicular to the current flow 
and lies in the (X, Y) plane (see Fig 1). The axis [Oil] 
and [233] define the two main crystallographic axis in the 
plane perpendicular to growth direction [311]. The am- 
plitude of each resonance exhibits a strong dependence 
on the intensity of B. In particular, with increasing B, 
the low- voltage resonances "A" decrease steadily in am- 
plitude, whereas the others, "B" and "C", have a non- 
monotonic magnetic field dependence. The Figure 3(b) 
and (c) show clearly two characteristic types of magnetic 
field dependence: type "A" shows a maximum on G(B) 
at B = T followed by a steady decay to zero at around 
8T; type "B" shows a broad maximum at ~ 4.5 T, fol- 
lowed by a gradual decay to zero. 
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FIG. 4. (a) I(V) characteristics in an in-plane magnetic 
field of 4.5 T. The first curve (circles) is for S|| [011]; the sec- 
ond curve (triangles) is for £>||[233]; (b) angular dependence 
of the peaks current. 

Figure 4(a) shows the I(V) characteristics in an in- 
plane magnetic field of 4.5 T. The first curve (circles) is 
for B 1 1 [011] ; the second curve (triangles) is for fl||[233]. 
We observe a clear anisotropy in the dependence of I(V) 
on B for the two field orientations. As can be clearly 
seen in Fig.4(a) peaks "A" and "B" in the I(V) plot 
reveal a strong anisotropy of about p ~ 0.5. We have 
also determined angular dependence of the peaks. The 
results are plotted in Fig. 4(b) for peaks "A" and "B". 
Note that all peaks, observed over the bias range (~200 
mV) have a maxima in current amplitude at orientation 
of a field B\\ [011] or B\\ [233]. The main effect to be noted 
from Fig. 4 is the dependence of the current as a function 
of the in-plane magnetic field orientation. 

We can understand the magnetic field dependence of 
the features in terms of the effect of B on a tunnelling 
electron. Let a, j3, and Z indicate, respectively, the di- 
rection of B, the direction normal to B in the growth 
plane (X, Y), and the normal to the tunnel barrier, re- 
spectively (see Figure 1 (b)). When an electron tunnels 
from the emitter into the dot, it acquires an additional 
in-plane momentum given by |l4| 



where As is the effective distance tunnelled along Z. 
This effect can be understood semiclassically in terms of 
the increased momentum along /?, which is acquired by 
the tunnelling electron due to the action of the Lorentz 



force. In terms of mapping out the spatial form of an 
electronic state, we can envisage the effect of this shift in 
as analogous to that of the displacement, in real space, of 
the atomic tip in a STM imaging measurement. The ap- 
plied voltage allows us to tune resonantly to the energy of 
a particular QD state. Then, by measuring the variation 
of the tunnel current with B, we can determine the size of 
the matrix element that governs the quantum transition 
of an electron as it tunnels from a state in the emitter 
layer into a QD. In our experiment, the tunnelling ma- 
trix element is most conveniently expressed in terms of 
the Fourier transforms of the conventional real 

space wavefunctions JT^Jl^p . Here the subscripts i and / 
indicate the initial (emitter) and final (QD) states of the 
tunnel transition. Relative to the strong spatial confine- 
ment in the QD, the initial emitter state has only weak 
spatial confinement. Hence, in fc-space corresponds to a 
sharply peaked function with a finite value only close to 
k = 0. Since the tunnel current is given by the square 
of the matrix element involving <&j(fc) and $QD(k), the 
narrow spread of k for $i (fc) allows us to determine the 
form of <&QD(k) by varying B and hence fc according to 
(|l|). Thus by plotting G(B) for a particular direction of 
B we can measure the dependence of \&QD(k)\ 2 along 
the fc-direction perpendicular to B. Then, by rotating 
B in the plane (X, Y) and making a series of measure- 
ments of 1(B) with B set at regular intervals ( A6 ~ 5°) 
of the rotation angle 8, we obtain a full spatial profile 
of \&QD(kx, ky)\ 2 ■ This represents the projection in fc- 
space of the probability density of a given electronic state 
confined in the QD. 

The model provides a simple explanation of the mag- 
netic field dependence of the resonant current features 
"A-C". In particular, the forbidden nature of the tun- 
nelling transition associated with "B" at B = T is due 
to the odd parity of the final state wavefunction, which 
corresponds to the first excited state of a QD. The ap- 
plied magnetic field (i.e. the Lorentz force) effectively 
breaks the mirror symmetry at B = and thus makes 
the transition allowed. 
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FIG. 5. Distribution in the plane (kx, ky) of the differen- 
tial conductance, G = dl/dV, for three representative states. 
This provides a spatial map of \&Qn(kx, ky)\ 2 , the square of 
the Fourier transform, &qd (kx, ky), of the probability den- 
sity of the electron confined in the dot. X and Y define the 
two main crystallographic axis, [Oil] and [233], respectively, 
in the (311)-oriented GaAs plane. 

Figure 5 shows the spatial form of G(B) ~ 
\&QD(kx,ky)\ 2 , in the plane (kx,ky) for the two rep- 
resentative QD states corresponding to the labelled fea- 
tures in Fig. 3(b) and (c). The fc- values are estimated 
from relation (1), assuming As has nominal value of 
30 nm which we estimate from capacitance measurements 
and from the doping profile and composition of the de- 
vice. The contour plots reveal clearly the characteristic 
form of the probability density distribution of a ground 
state orbital and the characteristic lobes of the higher en- 
ergy states of the QD. The electron wavefunctions have 
a biaxial symmetry in the growth plane, with axes cor- 
responding quite closely (within measurement error of 
15°) to the main crystallographic directions X — [Oil] 
and Y — [233] for (311) B-substrate orientation. For a 
similar In As QD structure grown on a (lOO)-substrate 
we also obtained characteristic probability density maps 
of ground and exited states. 

To summarise, we have observed features in /(V) cor- 
responding to resonant tunnelling through a limited num- 
ber of discrete states whose wavefunctions display the 
symmetry of the ground state and excited states of quan- 
tum dots. With the simple device configuration we have 
used, it is not possible to tell whether an excited state fea- 
ture and a ground state feature correspond to the same 
quantum dot. This question could be resolved by new 
experiments on structures with electrostatic gates. 

IV. CONCLUSION 

In conclusion, we have shown how the magnetotun- 
nelling spectroscopy provides a new means of probing the 
spatial form of the wavefunctions of electrons confined in 
quantum dots. The study revealed a biaxial symmetry 
of QD states in the growth plane. We observed the el- 
liptical shape of the ground state and the characteristic 
lobes of the higher energy states. 
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